Rigid Unit Modes in Tetrahedral Crystals 

Franz Wegner 

Institut fur Theoretische Physik, Universitat Heidelberg, 
Philosophenweg 19, D-69120 Heidelberg 

Abstract. The 'rigid unit mode' (RUM) model requires unit blocks, in our case 
tetrahedra of SiC>4 groups, to be rigid within first order of the displacements of 
the O-ions. The wave-vectors of the lattice vibrations, which obey this rigidity, are 
determined analytically. Lattices with inversion symmetry yield generically surfaces 
of RUMs in reciprocal space, whereas lattices without this symmetry yield generically 
lines of RUMs. Only in exceptional cases as in /3-quartz a surface of RUMs appears, if 
inversion symmetry is lacking. The occurence of planes and bending surfaces, straight 
and bent lines is discussed. Explicit calculations are performed for five modifications 
of SiC>2 crystals. 



PACS numbers: 63.20.-c 
1. Introduction 

Displacements of ions in a solid, which alter distances between neighboring ions, produce 
much stronger forces than those, which vary only angles between adjacing bonds. This 
has led to the idea of rigid-unit modes (RUMs), that is distortions, which do not change 
distances between the ions in a unit to first order in the displacements. Typical examples 
are crystals, in which silicon- or aluminium-ions are surrounded by four oxygen ions. 
The units of the tetrahedra of these oxygens are required to be rigid. Each oxygen ion 
belongs to two tetrahedra. 

There have been extensive numerical studies of RUMs of such crystals by Dove, 
Giddy, Hammonds and Heine et al. For a review see pQ. More recent presentations 
of RUMs in framework aluminosilicates can be found in the internet [2] and in the 
review[3]. These rigid- unit modes do not only signal soft phonon- modes, but they 
are also at the origin of a large number of displacive phase transitions. The first one 
considered was the transformation between a and (3 quartz [1]. Recent applications of 
RUM modelling include basic ideas for the development of new zeolites[5]. These zeolites 
are important catalysts in petrochemical refineries due to their high internal surface 
areas and molecular sieving properties. Another application deals with the flexibility of 
framework structures, which due to RUMs allows cation substitutions with a minimum 
of energy cost, since the geometric stress associated with the substitution is absorbed 
by rigid-unit type motion of the polyhedra near the substitution site [6]. 
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Basic to the numerical calculations is the computer program CRUSH [S] . In this 
program the rigid tetrahedra are assumed to be individual molecules, and harmonic 
forces are added between the two 'split' atoms, which should be one. It calculates the 
phonon frequencies o^q for given wave-vector q, which allows to determine a pseudo- 
intensity 

evaluated for small Q. For a wave-vector q with n RUMs this quantity approaches 
/(q) « n/Q. In this way the authors determined RUMs for a large number of crystals. 

There are also some analytic calculations along planes and lines of symmetry. Such 
a calculation has been performed by Vallade, Berge, and Dolino [9] for /3-quartz. The 
concept of RUMs itself goes back to Megaw [10] and Grimm and Dorner [I] . 

The present paper reports analytic calculations of RUMs. The variation of the 
lengths of the bonds of the units is calculated as a function of the displacements of 
the O-ions. Since the number of coordinates of the displacements is the same as the 
number of lengths of bonds, namely 6 times the number of tetrahedra, one has to find 
non-trivial solutions of linear homogeneous equations, which means the determination 
of the zeroes of the determinant of the corresponding coefficient matrix. In reciprocal 
space these reduce to 6n t equations, where n t is the number of tetrahedra in the unit 
cell of the crystal. So one has to calculate the determinant of a 6n t x 6n t matrix, where 
the wave vector q enters only through the three complex phase factors 

Pk = e ia * q , (2) 

where are the basis vectors of the lattice and q is the vector of the oscillating wave of 
displacements. We will show that for a crystal with inversion symmetry the determinant 
of an appropriately defined matrix is real for any given ps. Thus there is only one 
condition to be met with three unknowns, which in general defines surfaces of RUMs 
in reciprocal space. If inversion symmetry is lacking, then the determinant is complex, 
and two conditions have to be met for the RUMs: both the real and the imaginary part 
have to vanish. Correspondingly RUMs are found only at the intersections of the zeroes 
of the real and the imaginary part. Thus generically one obtains in these cases lines of 
RUMs, but not surfaces. It may happen, however, that both the real and the imaginary 
part have a factor in common, so that RUMs extend over a whole surface. An example 
is /3-quartz, where a whole plane of RUMs besides lines of RUMs are found. Although 
in many cases the RUMs are located on planes or lines, there are also cases where a 
surface can bend, as was already found for HP tridimyte by Dove et al [11]. Here (eq. 
l53l) the analytic expression for the surface is given. Lines of RUMs may lie in planes. 
But there are also non-planar lines. Such a line has been found in a-quartz (figure [6]). 

The general theory is given in section [2j Some remarks on planes and bending 
surfaces, straight and bent lines are given in subsection 12.11 The basic equations will be 
derived in subsection 12 .21 In particular the coefficient matrix for the determinant will 
be given. In subsection 12.31 the consequences of inversion symmetry will be considered. 
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It will be shown that the determinant multiplied by appropriate powers of pk is real for 
lattices with inversion symmetry. This section is concluded with a few remarks on the 
actual algebraic calculation in subsection 12.41 

This theory will be applied to five modifications of SiC>2 in section [3j first the 
crystals with inversion symmetry /3-cristobalite and HP tridymite are considered, then 
the RUMs for the crystals without inversion symmetry /3-quartz, a-cristobalite, and 
a-quartz are calculated. 



2. General idea 



2.1. Planes and bending surfaces, straight and bent lines 

The determinant (or its real and imaginary part) is factorized after calculation. The 
factorization procedure of algebraic computer programs is very useful, since the zeroes 
of the various factors yield different locations of RUMs. Often factors of the form 

=F AAA (3) 

with integer powers U and l\ appear. Then the determinant (or its real or imaginary 
part) vanishes along the plane 

[(h - ii)ax + (h - l' 2 )aa + (h ~ I'M ■ q = j ( (4) 

with integer m. If we introduce the reciprocal wave-vectors b defined by 

SLi ■ bj = 2tx 5 ij (5) 

and represent the wave-vector q in this basis, 

q = £ &bi = f bi + r]b 2 + Cb 3 , ai • q = 2vr^ (6) 

i 

then the plane is given by 

(h - i[)t + (h - + (is - QC = { ™ i (7) 

I m +2 

But there are also cases where the factor has a more complex form, which yields a 
winding surface or a winding line. Such a surface appears in HP tridymite, eq. (|53p . 
Lines of RUMs are planar, if they come from a factor of type of the real or imaginary 
part. If the factors of both the real and the imaginary part are of type ([3]) then the line 
of RUMs is straight. If none of it is of this type, then the lines are generically non-planar 
as in a-quartz (figure [BJ . 

2.2. Basic equations 

Suppose the oxygens are located at 

I"n,i = R-n,i + u n,i, (8) 

Rn,i = ^2(n a + c ia )a a = A(n + Ci) (9) 
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with the equilibrium positions R and the displacements u. The integers n a number 
the elementary cells. The atomic coordinates of the 2n t oxygens in the elementary cell 
are denoted by Ci i0l . The operator A maps the unit cube onto the elementary cell. It 
performs the similarity transformation 

Ax = ^a a (e Q • x), (10) 

a 

from the orthogonal unit vectors e a to the lattice vectors a Q . Thus the components of x 
and n + c expand the vectors in the basis {a}. The distance between the oxygens at the 
corners of the tetrahedra should be fixed in first order in the displacements. Thus for 
such a pair of atoms at r n+m j and r n+m /j with fixed distance one obtains the condition 

(R-n+m,i — R-n+m'j) ' (u n +m,i — u n+m' ,j) = 0. (11) 

The distance vector reads 

Rn+m,j Rn+m',j ^d m _ m ' j j , (1^) 

d ffl -m',i,j = Ci + m - m' - Cj. (13) 
Thus eq. (11 II) can be rewritten 

i^dm-m'iij ' (u n +m,i ^-n+m' ,j) 

dm— m' ' (A Un+m,i A Un+m'j) 
— dm-m'^j • (u n +m,i — Un+m',j) = 0, (14) 

where A T is the transposed operator of A, 

u := *4 T u = e a (a a • u). (15) 

a 

The operator A T allows to express 

^l T b a = 27ie a , A T q = 2tt£. (16) 
Since the vector d is independent of n the Fourier transform of equation (I14p yields 

d m -m',i,j • (pmUi(q) - Pm'Uj(q)) = (17) 



with 



k=l 



u «(q) = 5Z ex p( _i H n fc a fc ■ q)u n ,i- (19) 

n k 

The eqs. (TTTT) constitute a set of linear homogeneous equations for the displacements 
u. We are looking for those wavevectors q which allow nontrivial solutions of the set of 
equations i fTTj) . Note that A has disappeared from the eqs. (TTTT) . Thus lattice parameters 
of the unit cell do not enter into the calculation. Nevertheless inversion symmetry and 
rotational symmetries play still a role. To facilitate notation the 6n t edges of the n t 
tetrahedra are numbered by e. The edge e connects the vertices i e and j e . To this 
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edge the triples m e , m' e and the distance vector d e = d me _ m / j e j e are associated. Then 
equation (fTTl) can be written 

A e := d e • (p me ti ie (q) - p m iu je (q)) = 0. (20) 

This set of homogeneous equations in the elongations u has non-trivial solutions, if the 
determinant Ai of the matrix M with 6n t x 6n t elements 

M e>ka = — e = d eQ (p me 5j ejfc — p m 'Jj e ,kj • (21) 

vanishes. The index e of the edges denotes the rows, the indices k and a run 
independently from 1 to 2n t and from 1 to 3, resp., and number the columns of the 
matrix. 

In general the determinant of M yields a complex number as function of the p k . 
Thus both the real part and the imaginary part of M. have to vanish. The determinant 
can be expanded 

M := det(M) = ]T p tjkP \pip k 3 (22) 

ijk 

with real coefficients and integer i, j, k. Thus M.(p\, pi, pi) = Ai*(pi, P2, P3) holds. 
Therefore a vanishing M. for some wave- vector q also implies that it vanishes for — q. 

2.3. Inversion Symmetry 

We show that in the case of inversion symmetry M. multiplied by an appropriate factor 
p m is real. Suppose the center of inversion is located at 

R l = Ap = Y,P a * a (23) 

a 

and denote the sublattice obtained by inversion from the sublattice i by I(i). (i may, 
but need not be identical to I(i))- Then 

Cj + c m = 2p + \i (24) 

holds with integer If = If,*. Inversion is also applied to the edges, which generates /(e) 
from e. Then the corners at the ends of the edges obey 

Zj( e ) = I(i e ), jl(e) = I(je)- (25) 

Due to this inversion the two distance vectors add up to zero, 

d e + d /(e) = (26) 

with 

d e = c ie + m e - m' e - c je , (27) 
d/ (e ) = c/( ie) + m /(e) - m' /(e) - c /0e) . (28) 
Therefore one obtains 

s e := m e + mj( e ) + l ie = m' e + m' I{e) + l je , (29) 
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which yields 

m e + m/( e ) = (30) 

with 

m e = m e + lte - Se , rh /(e) = m /(e) + — - Se . (31) 
The matrix M is introduced by 

M e , ka = p Se /2P-l fc /2^fe,fca, (32) 
M e ^ a = d ea (pm e 5i e ,fc — Pm' e Sj e ,k) ■ (33) 

In rows e and /(e) the only non-zero matrix elements of M are 

■^e,i e a = <^eaPm e ; -^e,j e a = ^ecePm' e ^ (34) 

^/(e),i e c« = 4aP-m e = Ms*i e a; M I(e),j e a = ~deaP- m > £ = M e,j e a (35) 

Therefore exchanging all pairs of rows e and /(e) will transform the determinant M. 
into (— ) 3nt times the determinant in which all arguments pi are replaced by 1/pi (note 
that p_ m = 1/pm)- Since a tetrahedron does not transform into itself under inversion, 
inversion symmetry is only possible for even n%. Thus 

M{pi,p 2 ,p 3 ) =M[—,—,— ) ■ (36) 

\Pl P2 P3/ 

Since for real wave- vectors q one has 1/p, = p*, one deduces 

M = M*. (37) 

Thus M. is real for crystals with inversion symmetry. The connection between Ai and 
M. is obtained from eqs. (|3T| 1321 

:= det(M) = (psm)" 3/2 Al, (38) 

(p Sm )- 3/2 =n(psj- i/2 n(a) 3/2 - (39) 

e i 

The sum X) m has to be extended over all vertices of all tetrahedra. Since each vertex 
belongs to two tetrahedra, there are two contributions to each vertex. 

Although relation (|36l) holds for crystals with inversion center only, the 
transformation from M. to Af is also useful for lattices without inversion symmetry 
for the following reason: The introduction of the m is to some extend arbitrary. One 
may add an arbitrary m c to m at all corners of a given tetrahedron without changing the 
distance vectors d, m' = m + m c . Then in all six rows for the edges of this tetrahedron 
an extra factor p mc appears, M! = Pm c Af. Then one obtains X) m' = m + 4m c and 
Af = M.. Therefore Ai is invariant against this arbitrary choice in contrast to Al. The 
over-all sign of Al depends on the sequence of the vertices. Thus it is arbitrary. 

As a consequence crystals with inversion symmetry will show areas for RUMs in 
reciprocal space, since only one condition has to be fulfilled. If inversion symmetry is 
lacking, then both the real and the imaginary part of Al has to vanish. Unless both 
conditions coincide for whole areas, one obtains only lines for RUMs in reciprocal space. 
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2-4- Remarks on the calculation 

The calculation is performed by means of the algebraic computer program MAPLE. If 
the lattice is without inversion center, then M. is decomposed into its real 1Z\ and its 
imaginary part 7^-2- If besides the three factors pk there are no other variables, then 
the calculation of the determinant is extremely fast. If there is one extra variable like 
x in /3-quartz then it takes a few seconds. If there are three variables, like X2, y2, and 
Z2 in a-cristobalite and in a-quartz, then it takes to the order of hours. However if 
one assigns rational fractions like x<i = 23976/100000, then it runs very fast, whereas 
it takes quite a while, if one chooses decimal fractions like X2 = 0.23976. Apparently 
the exact calculation facilitates the calculation of the determinant, although it produces 
fractions with numerators and denominators of enormous size. Then also factorization 
still works, whereas it does not for decimal fractions due to rounding errors. 

The degeneracy of the modes can easily be calculated for RUMs on planes or 
on straight lines by determining the rank of the matrix M., since the corresponding 
restriction can easily be evaluated. If surfaces or lines are bent, then this is more 
difficult, since the constrained pu have to be given explicetely. In the following we will 
not determine the degeneracies. 

The crystallographic data used here are those of [12] . They describe average 
positions of the ions. They are almost not representative of the real local structure. 
Indeed the average positions are inconsistent with pair distribution functions [13]. This 
obviously is due to the RUMs. Thus the present analysis gives the RUM spectrum of 
the average structure. 

To the extend the location of the RUMs is given by symmetry, there will be obvious 
agreement between those determined in refs. [H EJ [3] and ours. For bent lines many 
agreements will be found. In comparison three things have to be kept in mind: First 
of all the cited references give only graphical representations. Therefore comparison is 
made by appearance only. Secondly, unless the atomic coordinates are determined by 
symmetry, theirs may differ from the ones used here. Thirdly, apparently Dove et al have 
included quasi RUMs that is regions in reciprocal space with a low phonon frequency, 
whereas here only zero frequency modes are taken into account. 

3. Various SiC>2 Crystals 

In the following the RUMs for five modifications of SiC>2 crystals are determined. 
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The coordinates c' for the Si-ions and c for the O-ions are given by 









43 


i 




Cj,2 


c i,3 


i 




1 











1 


1/8 


1/8 


1/8 


2 


1/2 








2 


-1/8 


-1/8 


-1/8 


3 





1/2 













4 








1/2 



ai = f e 2 + f e 3 
a 2 = f e x + f e 3 
a 3 = fei + §e 2 



The corners of the two tetrahedra are obtained by applying A to 

1st tetrahedron: c 1; c 2 , c 3 , c 4 ; 

2nd tetrahedron: ci, c 2 — ei, c 3 — e 2 , c 4 — e 3 . 

The determinant evaluates to 

(1 - Pl )(l - p 2 )(l - p 3 )( Pl - p 2 )( Pl - p 3 )(p 2 - p 3 ) 



(40) 



With eq. 



M 

one obtains 
1 ~ Pk 

1/2 
Pk 

Pk - Pi 



3/2 3/2 3/2 
2 12 Pl P2 P 3 

- 2isin(a fc q/2) = -2i sin(7r£ fc ), 
2isin((a fc - a z )q/2) = 2isin(7r(^ - &))• 



(PkPi) 1 / 2 

Thus all RUMs are located in the planes in reciprocal space 

(0,77,0, (£,0,0, (e^o), (e,e,o, (£,?7,o, (^^^)- 

The crystal has a center of inversion. One easily checks eq. ( 136|) for this crystal. 



(41) 

(42) 
(43) 

(44) 



Comparison The planes of RUMs found agree completely with those given in ref. pQ. 
3.2. HP tridymite 

The atomic coordinates of the HP tridymite phase of Si0 2 are given by 
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a 2 = |aei + ^ae 2 
a 3 = ce 3 



(45) 



The eqs. ffTTI) have to be fulfilled for the edges of the four tetrahedra with corners 

obtained by applying A to 

1st tetrahedron: c 1; c 3 + e 2 , c 4 , c 5 ; 
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2nd tetrahedron: c 2 — e 3 , c 3 , c 4 + e 1; c 5 ; 
3rd tetrahedron: c 2 , c 6 , c 7 + e 1; c 8 ; 
4th tetrahedron: ci, c 6 + e 2 , c 7 , eg. 
The determinant Ai yields 



n = 9 (p 3 ±£) -4/ 1 (p 1 ,p 2 ) + 14, (47) 

;i + p 1 )(i + p 2 )(i + p lP2 ) 



h{px,P2) 



/2(P1,P 2 ) 



P1P2 

8 cos(axq/2) cos(a 2 q/2) cos((a! + a 2 )q/2) 
= 2(1 + cos(aiq) + cos(a 2 q) + cos((ai + a 2 )q)), (48) 
(1 -pi)(l -p 2 )(l -P1P2) 



P1P2 

= 8isin(aiq/2) sin(a 2 q/2) sin((ai + a 2 )q/2) 

= 2i( sin(aiq) + sin(a 2 q) - sin((ai + a 2 )q)). (49) 

These functions fi and / 2 as well as the later introduced functions fs, eq. (1651) . and f^, 
eq. (11081) are invariant under rotations by 27r/3 around the z-axis, which transforms 

a i — > a 2 — > — a i — a 2 — > a i, 

Pi -> P2 -> 1/ (P1P2) -> Pi- (50) 
Due to the factor 1 — p 3 in eq. (j4"S|) RUMs are located in the plane 

(£,?7,0). (51) 
The factors of / 2 yield RUMs in the planes 

(0,7^,0, (£,0,C), (£,-£,0- (52) 

In contrast the zeroes of 7£ describe a winding surface in reciprocal space, which may 
be written 

16 7 
cos(27r£) = — cos(7r£) cos(7rr/) cos(7r(£ + 77)) — - 
9 9 

4 1 
= - (cos(2tt0 + cos(2vrr7) + cos(2vr(£ + 77))) - -. (53) 
9 3 

Note that the maximum of the sum of the r.h.s is +1, which is obtained at [0,0,0]. The 

minimum of the r.h.s. is —1, which is reached at [±1/3, ±1/3, 1/2]. For £ = 1/2 or 

77 = 1/2 or £ + 77 = 1/2 one obtains cos(2ttC) = -7/9, which yields ( = 1 = ±0.39183. 



Comparison These results are in full agreement with those obtained by Dove et al. 
[TT] by means of their numerical CRUSH-program [3, E] and reported in ref. pQ. 
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Other derivation One may determine this bending RUM also in the following way: One 
starts with the equations for the edges between c 3 , c 4 , and c 5 . Since they all lie in the 
xy-plane, the third component does not enter and one obtains the equations 



Pi«4,l = 1*5,1) 
p 2 W 3j2 = #5,2, 

#3,1 + #3,2 = Pi (#4,1 + W4,2) • 



U4,l = Us,l, 
#3,2 = #5,2, 

P2(£t 3 ,l + #3,2) = #4,1 + #4,2, #3,1 + #3,2 = Pi (#4,1 + #4,2)- (54) 

Evidently, if 

Pi ^ 1, P2 ^ !, P1P2 ^ !, (55) 
then all these components vanish 

#3,1 = #3,2 = "4,1 = #4,2 = #5,1 = #5,2 = 0. (56) 

Similarly one shows by considering the third and fourth tetrahedron that under the 
same condition (1551) one obtains 

#6,1 = #6,2 = #7,1 = #7,2 = W 8 ,l = #8,2 = 0. (57) 

This implies that the tetrahedra are rotated around axes parallel to the xy-plane. There 
are twelve equations left for the twelve other components u. One can use six of them to 
eliminate w 3j3 to Us,3- Finally one calculates the determinant of the coefficient matrix 
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(5f 



of the last six equations for the six components of u x and u 2 , which has to vanish. This 
determinant factorizes in a factor p 3 — 1 and into 




P1P2 



9/1 



= 4( cos(a!q) + cos(a 2 q) + cos((ax + a 2 )q)) — 9 cos(a 3 q) — 3 
which again yields eq. floH . 



0, (59) 



3.3. (3 -Quartz 

The coordinates for the Si and O-ions are given by 
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(60) 
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The corners of the three tetrahedra are obtained by applying A to 

1st tetrahedron: c 2 + e 1; c 3 — e 3 , c 5 , c 6 + ei — e 3 ; 

2nd tetrahedron: ci, c 3 + e 2 , c 4 + e 2 , c 6 ; 

3rd tetrahedron: c 1; c 2 + e x + e 2 , c 4 + e x + e 2 , c 5 . 

The determinant reads 

m = [—^ — -J -jj^{n x +K 2 ), (ei) 
ftl = ^-^7^-/2(^1,^2) 

Ps 

x (-k!(^^)+k l + k 2 f 1 (p 1 ,p 2 )), (62) 

1 — Di 

TZ 2 =2x(2x-l)(3x-l)(4x-l)(6x-l)^ 7 p/ 3 (p 1 ,p 2 ), (63) 

Ps 

Mp^)= (fil - P2)(1 -?f {1 -" lP2 \ (64) 

PlP2 

fci = 2x(4x- l)(432x 4 - 540x 3 + 252x 2 - 51x + 4), (65) 

k 2 = -x(4x- l)(3x- l) 2 , (66) 

A; 3 = 12x 2 - 6x + 1. (67) 

Due to the factor (1 — p 3 ) RUMs are obtained in the plane 

(£,?7,0). (68) 

The factor (1 + p 3 ) in and the zeroes of / 3 in 7£ 2 yield RUMs along the straight lines 

[£,£,1/2], [£,-2£,l/2], [-2£, £,1/2]. (69) 

The zeroes of f 2 in IZi and of / 3 in 7£ 2 yield RUMs along the straight lines 

[0,0, C], [0,1/2,C], [1/2, 0,C], [1/2,1/2,0- (70) 

Finally the zeroes of the large parenthesis in the expression for IZi and the zeroes of / 3 
in 1Z 2 yield RUMs along the curves 

[£,£,C(0], [£,-2£,C(0], [-2£,£,C(£)], (71) 

where £(£) is given by 

ki + 8k 2 cos 2 (7r£) cos(27r£) 

>3 

With x = 0.4202 of ref. [12] one obtains 

fei = 0.26799, fc 2 = -0.019427, k 3 = 0.59762 (73) 

and 

cos(2ttC(£)) = °- 6278 - °- 3641 cos 2 (vrO cos(2tt£), (74) 



cos(2ttC(0) = ^3 — ■ (72) 
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However, this value of x yields tetrahedra, which are far from being equilateral. The 
value for equilateral tetrahedra is 
1 1 



x 



0.2113, 



2 V12 
which yields 

k x = -0.015801, k 2 = 0.004380, k 3 = 0.26795, 
cos(2ttC(0) = -0.41068 + 0.91068 cos 2 «) cos(2<). 
The function is plotted for both values x in fig. [U 



l /2 



x=0.4202 



+ '/2 



+ 72 -72 



l /2 



x=0.2113 



+ 72 



'/2 



+ 72 



(75) 

(76) 
(77) 



Figure 1. of the RUMs of /3-quartz given by eqs. fTT j I72 |) for x = 0.4202 and 
x = 0.2113. 



Comparison We find agreement of the RUMs for /3-quartz in reciprocal space with 
those given in refs. [TTJ [H [21 [3]. This includes the winding lines of fig. [1] seen in ref. 
[2, 3J with the ideal value x of eq. ( 1751) . 

3.4- a-Cristobalite 



The coordinates d for the Si-ions and c for the O-ions are given by 











i 


Ci,l 


Ci,2 


Ci,3 














1 


X2 


2/2 


^2 






% 


<i 






2 


-x 2 


-2/2 


1/2 + z 2 






1 









3 


1/2 - y 2 


1/2 + x 2 


1/4 + 2 2 


ai 


= aei 


2 


— Xi 


— Xi 


1/2 


4 


1/2 + 2/2 


1/2 - x 2 


3/4 + z 2 


a 2 


= ae 2 


3 


1/2 - xx 


1/2 + x x 


1/4 


5 


2/2 


x 2 


-z 2 


a 3 


= ce 3 


4 


1/2 + 3* 


1/2 - xi 


3/4 


6 


-2/2 


—x 2 


1/2 - z 2 














7 


1/2 - x 2 


1/2 + 2/2 


1/4 - z 2 














8 


1/2 + x 2 


1/2 - 2/2 


3/4 - z 2 







The corners of the four tetrahedra are obtained by applying A to 

1st tetrahedron: c 1; c 4 — e 3 , c 5 , c 7 ; 

2nd tetrahedron: c 2 , c 3 — ei — e 2 , c 6 , c 8 — ei — e 2 ; 
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3rd tetrahedron: c x + e 2 , c 3 , c 6 + e 2 , c 7 ; 
4th tetrahedron: c 2 + e 1; c 4 , c 5 + ej + e 3 , c 8 . 
The determinant reads 



^2 



fc 4 (^l+^ 2 ), 
2 

5^5, cos'(2vre 3 ), 

„ (1-P1)(1-^)(1-P|)(P1-P2)(1-P1P2) 

9 9 

P1P2P3 

^ Afyi cos i (2< 1 ) cos , '(27r6) J 



A; = -(2x 2 - l)(8z 2 y 2 + x 2 - y 2 ) 



with 



fc 



(79) 
(80) 

(81) 
(82) 

(83) 
(84) 



The sum over i and j runs at most up to i = 3, j = 3, i + j = 4. The coefficients k^i 
are polynomials in x 2 , y 2 , z 2 . They are of order 4 in z 2 and run up to order 8 in x 2 and 
y 2 . Thus in general they are lengthy expressions. 
For 



x 2 = 0.23976, y 2 = 0.10324, z 2 = 0.17844 
from ref. [12] one obtains 

fc' = -0.00052747, k = -0.018470, 
and the non-zero cofficients k^i 



(85) 
(86) 



% 


j 


lOO^o 








-8.0953 





1 


9.6093 





2 


9.4540 





3 


-10.9503 


1 


1 


-8.4081 


1 


2 


-7.3564 


1 


3 


7.9854 



i 


j 


100%! 








-0.7104 





1 


-3.7422 





2 


2.5622 


1 


1 


-19.5336 


1 


2 


18.3067 


2 


2 


-15.9707 



i 


j 


100fc ii2 








8.8057 





1 


-5.8672 


1 


1 


3.9092 



(87) 



The zeroes p\ = ±1, p 2 = ±1, p 2 = 
of TZi from 

-(ft ± \fd 



Pi, p 2 = 1/pi of 7?. 2 yield RUMs for the zeroes 



cos(2ttC(0) 



2^2 



d 



9\ 



40002- 



For pi = 1 and p 2 = 1 one obtains RUMs along lines 

[o,e,c(0], [e,o,c(0], 



(88) 
(89) 
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resp. with 

10 2 g 2 = 2.9385 - 1.9579 cos(27r£), 

10 2 #i = - 1.8904 - 4.9691 cos(2ttO + 4.8982 cos 2 (27rO, 
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10 4 d = (1 - 
They are shown in fig. [21 



cos(2tt£)) 2 (3.3647 + 4.1431 cos(2tt£) + 0.7718 cos 2 (2vr£)).(90) 



Cf+V2 




+ V2 



Figure 2. of the RUMs of a-cristobalite described by eqs. (j88l l89 l [90]) . 



For pi 



resp. with 



— 1 and p 2 = — 1 one obtains RUMs along lines 

[i/2,£,c(£)], K,iAC(fl] 

10 2 # 2 = 14.6729 - 9.7764 cos(2tt£), 



10V = 5.5940 + 34.0981 cos(2tt£) - 31.7152 cos 2 (2ttO, 
10 4 c/ = - 127.156 + 330.023 cos(2<) - 74.142 cos 2 (2vr£) 
- 394.118 cos 3 (2ttO + 265.367 cos 4 (2<). 

The calculation with (|85|) did not give real solutions C(£)- 
For p 2 = pi and p 2 = 1/pi one obtains RUMs along lines 

[£,e,c(0], [£,-£, c(£)] 

with 



(91) 



(92) 



(93) 



10 2 # 2 = (2.9674- 1.9772 cos(20) 2 , 

lO 2 ^ = - 0.7104 - 7.4844 cos(2<) - 14.4092 cos 2 (2^) 

+ 36.6135 cos 3 (2tt£) - 15.9707 cos 4 (2ttO, 
d =(2g 2 + 9l ) 2 . (94) 

The square root of d is rational in this case. The first solution (1881) yields the straight 
line of RUMs 



[£,±£,o], 



(95) 
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The second one yields the relation 
cos(2ttC(0) : 



92 



15 



(96) 



The locations of these RUMs in reciprocal space are shown in fig. [3J 



Cf+'/2 




+ 72 



Figure 3. of the RUMs of a-cristobalite given by eqs. ([Ml I? 

The RUMs given by the zero p 3 = 1 are the RUMs eq. 0951) just mentioned plus 
RUMs along the lines [£, rj, 0] with £ and 77 related by 

cos(2vr£) cos(27r7?) - A(cos(2tt£) + cos^vrr?)) + B = 0, (97) 
where for the above given £2, 1/2, ^2 one obtains 

A = 1.3713, B = 1.5048. (98) 
The corresponding line is plotted in fig. HI 




+ 72 



Figure 4. 77 versus £ of the RUMs of a-cristobalite given by eqs. (|97l |98|) . 

Finally the RUMs given by the zero p 3 = —1 of 7I2 are given by lines [£,77, 1/2] 
with £ and rj related by 



^ fc i jCOS i (27r£) cos J '(2vrr/) = 0, 

M=0 



(99) 
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which with (I85p yields 

15.9707 cos 2 (2ttO cos 2 (2vrr/) 

+ 7.9854 cos(2ttO cos(27n?)(cos 2 (20 + cos 2 (27rr/)) 

— 25.6631 cos(27r£) cos(27r?7)(cos(27r£) + cos(27r^)) 

- 10.9503(cos 3 (2ttO + cos 3 (27r7?)) 

+ 15.0347 cos(2<) cos(27r??) + 6.8918(cos 2 (2<) + cos 2 (27rr/)) 

+ 7.4844(cos(2<) + cos^vrr?)) + 1.4208 = 0. (100) 



This can be rewritten 



4 - 4.5851c 3 + 3.6089c 2 + 3.7491c s + 0.7117 
c 2 + 0.9001c s + 0.1567 



(101) 



with 



C s = COs(27r£) + COs(27T77), Cd = cos(27r£) — cos(27r?7). (102) 

The numerical calculation showed that either c\ from eq. fllOlH is negative or one of the 
cosines is larger than one. Therefore there is no real solution of this type. 



Comparison The figure for a-cristobalite in ref. [21 E] shows the RUMs of fig. [2] and 
[Hand the line eq. [95] (green) for £,?7,C > 0. The line fig. [3] cannot be seen. However, 
a region of quasi-RUMs (modes of small frequency) is shown, which probably hides the 
line of fig. [31 



3.5. a -Quartz 

The coordinates of the Si- and O-ions in a-quartz are 



% 


4i 


<2 


<3 


1 


Xi 





2/3 


2 







1/3 


3 


-Xi 


—X\ 






i 


Q,i 


Cj,2 


Ci,3 






1 


X 2 


2/2 


z 2 






2 


~V2 


x 2 - 2/2 


2/3 + z 2 


ai 


= \ae l - 


3 


2/2 - x 2 


— x 2 


1/3 + z 2 


a 2 


= \ae l + 


4 


2/2 


x 2 


-z 2 


a 3 


= ce 3 


5 


-%2 


2/2 - x 2 


2/3 - z 2 






6 


x 2 ~ 2/2 


-2/2 


1/3 - z 2 







2 C 



The corners of the three tetrahedra are obtained by applying A to 
1st tetrahedron: ci, c 2 + e x — e 3 , c 5 + e x + e 3 , c 6 + e 3 ; 
2nd tetrahedron: c 2 - e 3 , c 3 + e 2 - e 3 , c 4 + e 3 , c 6 + e 2 + e 3 ; 
3rd tetrahedron: Ci — e 3 , c 3 + e 1 + e 2 — e 3 , c 4 + e 3 , c 5 + e x + e 2 . 
The determinant evaluates to 

= ^(7^ + 7*2), 

2(3a; 2 - 2)(— 2x 2 + 3z 2 x 2 - 6y 2 z 2 + 5y 2 ) 



M 

h 



/2(Pl,P 2 )" 



Pl 



Pa 



9 

k 2 fi{pi,p 2 ) + k 



1 + 

P3 



+ fc 4 



(104) 
(105) 

(106) 
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-2 — (.fafl(j>l,P2) + h) 
P3 



+ 



i-pl 

P3 



(k 7 fi(pup2) + k 8 fi(pl,pl) + k 9 f 4 (p u p 2 ) + ho) 



1 + p 2 

+ (hi — — + k 12 )f 3 { Pl ,p 2 ), 

P3 

f( (Pl + P2)(l + P1P 2 2 )(1+ P{P2) 
H\Pl,P2) — o~2 • 



P\p\ 



(107) 

(108) 



with lengthy polynomial expressions of x 2 , y 2 , z 2 for the constants k 2 to kx 2 . With 

x 2 = 0.4141, y 2 = 0.2681, z 2 = 0.7854 (109) 
from ref. [12] one obtains 



(110) 



jfej = -0.037819, k 2 = 0.018997, k 3 = -0.030172, 

k 4 = 0.000088, k 5 = 0.006248, k 6 = -0.060332, 

k 7 = 0.021768, k 8 = -0.009082, k 9 = -0.000285, 

k 10 = -0.078514, k u = -0.006089, k 12 = 0.001074. 

From the zeroes of TZ\ combined with the zeroes of TZ 2 one obtains three classes of 
solutions: 

(i) f2(pii P2) — yields a class of RUMs 

[£,o,c(OL [o,-£,C(OL K,£,C(0] 

where is determined by 



(111) 



MPi 



p 2 = l,p 3 = e 



0. 



(112) 



The solutions are shown in fig. [5j One line of RUMs is nearly straight, the other one is 
snaky. 



St + 72 




+ 72 



-72 



Figure 5. of the RUMs of a-quartz described by eqs. (|lllQ112p . 
(ii) 1 — p\ = and f$(pi, p 2 ) = yield straight lines of RUMs 

[e,e,o], [e,-2e,o], [-2^,0], 
[e,e,i/2], [e,-2e,i/2], [-2c e, 1/2]. 



(113) 
(114) 
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(iii) The condition 

k 2 fi(pi,P2) + h(Ps + i/p 3 ) + h = o, n 2 = o (lis) 

yields the curve depicted in fig. EJ This curve is not a planar curve. The projection 
onto the £,77) plane is nearly a circle. In £ direction it oscillates between —0.0433 and 
+0.0433. Thus it stays close to the ( = plane. Approximately it can be described as 
function of the parameter x 

7T 7T 

£ = r(x)sin(- - x), V = Kx) sia (- + X), (H6) 

o 

r(x) = 0.3290 - 0.0035 cos(6x), C = 0.0452 sin(3x) + 0.0019 sin(9x). 

This curve intersects with the lines of RUMs f 1 1 1 3 j) at \ = kn/3, and with the RUMs 
given by the curves (11111) at x — 71 / 6 + kn/3 for integer k. 




Figure 6. The RUMs given by eq. (|1 15[) shown in the £, 77-plane and C versus £ + r\ 
and £ — 77, resp., for a-quartz. 



Comparison Comparing the distribution of the RUM- vectors for a-quartz given in [21 [3] 
with the curves calculated here, we find agreement for the straight lines and the curves 
of figure dSJ). We do not find the non-planar curve figure El eq. (11161) . Since fig. 10 
of [3] shows only RUMs for positive (, this line should run through the dark regions of 
quasi-RUMs in the figure. 
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